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The reflection coefficient for the basic mode in a widening, straight, two-dimensional 
waveguide is computed for small wave numbers by using the perturbation method with 
the electrostatic case as the unperturbed case. The problem is treated as a perturbed infinite 
system of inhomogeneous linear equations, and it is shown that the matrix of the unperturbed 
system (which corresponds to the electrostatic case) can be inverted explicitly by using 
conformal mappings and physically unrealistic modes. Questions of convergence are 
discussed, and other examples for application of the method are indicated. 



1 . The Physical Problem 

1.1. Structure 

The problem which we are considering is a two- 
dimensional one arising from a three-dimensional 
waveguide structure in which the perfectly conduct- 
ing surfaces extend from — oo to °° in the direction 
of the s-axis in a Cartesian coordinate system. 
The intersection of these conducting surfaces with 
the #,?/-plane is given by the six lines: 

x>0, y = qir (0<g<l) 

x>0, y=—qir 
x=0, qj<y<ir 
x=0, —qir>y>—T 
x<0, y=ir 
x<0 y y=—Tr. 

The segment of the y-axis between — qw and qw, 
which we shall refer to as the aperture, separates 
the waveguide structure into two simple regions, 
I and II. 

1.2. Conditions on the Electric Field 

Let the electric field E=(E x ,E y , E z ), and assume 
that the time dependency is given by e iut where 

o)=kc, k=wave number, c= velocity of light. 
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Then we have the conditions 

(1) E x , E y , and E z satisfy Au + k 2 u = 0; 

^ as a£ 

1 ; dx + dy + dz ' 

(3) the tangential component of E vanishes at 

the boundaries. 

We shall also make the following assumptions: 

(4) E y is an even function of ?/ 

(5) E z =0. 

-> 

The behavior of E at infinity must satisfy 

E x -*0 



(6) lim \E y -TE e Ux \=Q, 



(7) lim \E y —E o e m —pE o e- ikx \=0, 

X^>+ a> 



E x -M) 



(8) I j {E 2 x -\-E 2 y )dxdy is bounded over every 
finite region. 

We shall introduce the following definitions: 

r (condition (6)) is called the transmission 
coefficient. 

p (condition (7)) is called the reflection coeffi- 
cient. 

E e ikx is called the incident wave. 

E y —E n e ikx is called the diffracted wave. 
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We shall now show that the field is completely 
determined by E y alone. Conditions (2) and (5) 
give 



(20 



bEy 

by 



*E X 



This shows that E y determines E x up to an arbitrary 
function of x only. But since E x and E y satisfy the 
reduced wave equation Au J rk 2 u=0, the same must 
be true for both sides of (2'). Furthermore, from 
(6) and (7) we know that E x must approach at oo . 
This fact gives us the uniqueness we seek. For if 
(2'), with a fixed E y , gave rise to two different 
solutions for E x , then their difference could only be 
a function of x alone. If this difference is D(x), 
then 



dx 2 



D(x)+k 2 D(x)=0 



and 



Z>(x)->0 at oo. 



This implies D(x) = 0. Thus, in what follows, we 
shall limit our attention to the determination of the 
^/-component of the electric field. Our first task is 
to write the bound ar}^ conditions on E y . 

From (1) we have 

(a) AE y 4-k 2 E y =0. 
From (3) we have 

(b) ^(x,±qw)=0 x>0 

(c) ^fe±T)=0 x<0 



(d) E y (o 7 y) = 



q.Tr<\y\<Tr. 



Note that in (b) and (c) we used the fact that E x =0 
everywhere on those boundaries. This implied that 

|X Tp 

-^=0, there; and then (2') was used to get (b) and 
(c). Continuity conditions in the aperture are 

(e) E y (0+,y)=E y (0- } y) \y\<qT 

(f) ^(0 + ,»)=^(0- f ») \y\<qr. 

And from (6) and (7) 

(g) Wm\E y -rE o e m \=0 

(h) lim \E v -E e ikx -pE o e- ikx \=Q. 



2. Method and Summary 

The problem described above has been treated be- 
fore, and the material has been presented by Marcu- 
vitz [1951, p. 141] and by Saxon [1943]. Quantita- 
tive data have been given by Marcuvitz [1951, p. 
307]. Furthermore, the problem has been treated 
earlier as a perturbation of the electrostatic case 
[Marcuvitz, 1951, p. 153], a method which will be 
used as the starting point for the present paper also. 

The essential feature of the present approach is 
the explicit inversion of an infinite matrix which 
characterizes the electrostatic case. In the special 
case where Q—% } this has already been done by 
Magnus and Oberhettinger [1950], who used alge- 
braic relations connecting the matrix elements for 
this purpose. The present paper uses instead, inte- 
gration in the complex plane, an approach which is 
of much wider applicability. The details may be 
described as follows: 

The first step toward finding the reflection coeffi- 
cient at the interface is to expand the field in each 
of the regions (I and II) in Fourier series. The modes 
are determined on each side by applying all the 
boundary conditions except the matching condition 
at the interface. The solution is then assumed to be 
a series in these modes with constant coefficients. 

Applying the matching conditions at the interface 
gives rise to two infinite sets of linear equations, one 
arising from matching the fields, and the other from 
matching the normal derivatives. Using these 
two sets of equations we can eliminate one group of 
coefficients, leaving an infinite matrix equation for 
the other group and a scalar side condition involving 
the reflection coefficient. Solving this matrix equa- 
tion, then, will allow us to apply the scalar side con- 
dition, and this will yield the value of the reflection 
coefficient. 

The matrix equation which arises from the match- 
ing conditions is derived in full and is valid for all 
values of the wave number, &<1. The coefficients 
of the expansion of the field in region II give the 
unknown vector. The equation is solved using a 
perturbation method. The electrostatic case (k = 0) 
serves as the unperturbed case. We shall see that 
the ability to solve the electrostatic case for arbi- 
trary right-hand sides provides us with enough power 
to determine the higher order terms in the electro- 
magnetic case (Jc5*0). Using this fact, and the scalar 
side condition, we shall determine the reflection co- 
efficient up to terms of order k 2 . 

The one remaining facet of the method is the proof 
of solvability of the electrostatic case for arbitrary 
right-hand sides. Since the equation in this case is 
Laplace's equation, we may use conformal mapping 
as an aid. We map the whole waveguide onto the 
infinite strip, using the Schwartz-Chris toff el formula. 
This enables us to make use of the fact that each un- 
known is a Fourier coefficient which can be written 
as an integral (across the aperture) of the field. The 
field is written as the derivative of a harmonic po- 
tential. In this way the matrix equation may be 
solved in the electrostatic case for a specific right- 
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hand side. However, this same technique can be 
used for more general right-hand sides. In order to 
accomplish this, we introduce into the original prob- 
lem "unrealistic" modes, i.e., modes which do not 
die out at +00. The iVth unrealistic mode will give 
rise to a right-hand side consisting of JV+1 nonzero 
entries. Thus, solving for all such unrealistic fields 
will yield a solution matrix whose product with the 
original matrix is a triangular matrix. The entries 
in the triangular matrix can be computed, using the 
conformal mapping; thus, the original matrix can be 
inverted. Therefore, we shall have solved the elec- 
trostatic case for arbitrary right-hand sides. 

3. Infinite System of Linear Equations 

In all that follows we shall use the notation that 
the restriction of E y to region I is denoted by E\ and 
similarly for E y l . When referring to the entire wave- 
guide we shall write simply E y . Inasmuch as our 
ultimate goal is to apply a perturbation technique to 
this problem, we shall assume that &<1 so that 

n 2 — F^ 2 >0 for all integers n^O. 

By separating variables in (a) and applying condi- 
tions (b) and (h) in region I we have 



k n 
i—x 



El(x,y)=:E [e ikx +pe-^}+j: c n e l T x cos -y (1) 



where 



k n =i-\n 2 - 



.k 2 gf, ti-1,2,.. 



In region II we have the additional condition (d) 
placed on the modes. So we must introduce a set of 
functions, <j> N , satisfying (a), (c), (d), and (g). In 
addition, let the JVth function satisfy 

N 
(e') ct> N (Q-,y)=cos~y, \y\<qTr 

in the aperture. The reason for condition (e') is 
that the wave coming in from the right is a super- 
position of cosines in the aperture. 

It will be necessary for us to have the 4> N explicitly, 
so we shall derive them here. The JVth mode must 
satisfy the conditions indicated in the sketch below. 

Separating variables, we see that <t> N can be ex- 
panded in a series of the form 



4>n(x, 2/)=S a ( n Y) e ilnX cos ny 

n=0 



1 ecu ion II) 



where 



r 



i-\!n 2 —k 2 n?*0 
k n=Q. 



N 1 
cos my cos — ydy. 



At x=0 we have 

*»(o, If)-ij«.« cos ny=\^T 1 y M<«* 

lo jpr<MO. 

Consider the case iV^O, and multiply both sides by 
cos my and integrate 

1 /•«» N 
C=- cos — ydy=0 

IT Jo (1 ' ' 

9 C q 

7T Jo 

Thus 
4to(fi, 2/)=S S mN e il ^ cos my N=l, 2, . . . (2) 

771 = 1 

where, for m= 1,2, 3, . . . and N=0, 1, 2, . . . 

9 r<jTr at 

Similarly, for N=0, we get 

<t>o(x, y) = qe ikx +^2 S mo e il m x cos my. (3) 

ra = l 

The field in region II may now be written 

E?(x,y)=ihd n <l> n (x,y). (4) 

n=0 

Applying condition (e), we have 



2 f°* N j , lW 2msin (qmr) 

cos my cos — ydy=(—\r —, — «, — rfcy-^r 

it J J (i JJ ; Tr(m 2 —N 2 /(i 2 ) 



# (l+p)+S Cn COS - 2/ = S dm<i>m(0, y) 
n = l <£ m=0 



\y\<<pr 



but, b} r definition, 



m 



<t> m (o,y) = cos-i 



\y\< ( iTT' 



Therefore, 

E o (i+ P )=d 

c n =d n n=l,2, 
Differentiate (1) with respect to x 
*El(x,y) 



(5) 



dx 



Let 



ikE [e ikx —pe ikx \+-y^c n k n e q cos-y. 
q_n=\ <r 

&l>m(x, y) 



$m(x,y) = 



dx 



Condition (f) gives 
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1 - °°. 77, 

Solving for the Fourier coefficients on the left-hand 
side, 

1 fQr oo 

ikE {\-p)=— J±d m + m {o,y)dy (6) 



2 r Qv m u 

Cn=-j77 H d m \l/ m (o,y) cos -ydy 



w=l f 2, 



(7) 



Combining (5) with (7), 



</ 



* _ dfc. 



> /^27T CO 

if S « 



_d m \l/ m (o,y) cos -ydy 
=o ff 



»=1,2, 



From (2) and (3) it can be seen that 

to(%,y) = ik(ie iJcx +i S S T0 l T e il r x cos n/ 

r = l 

tm(x,y) = ilt, l r S rm e il r x cos ry m=l,2, . 

r = l 

Therefore, 



(8) 



(9a) 



(9b) 



S (W»(c, y)=i&2d,+i XI <4 S ^m^r cos ry. 

m=0 m=0 r=l 



(10) 



Substituting into (8), 



kqd + ^i d m S ^SVm COS H/ I COS - ^ 
ra=0 r=l J ^ 



<}• 



»=1.2, 



co co 7 

to=0 r=l iCn 



71=1,2,.. 



(ii) 



Equation (11) is an infinite system of linear equations 
in the unknown Fourier coefficients, d n . 

A scalar side condition is obtained by substituting 
(10) into (6). 

I f*QTT C CO CO *~\ 

ikE {\ — p) =— 4 ikqd -\-i ]C d m ]C # m Z r cos ry V dy 
QTrJo L w=o r=i J 

/£ oo oo ^ 

=ikqd -\— X) dm 2] - # m Zr sin r7rg. (12) 

qw m =Q r =l f 

It can be seen from (3), (5), (4) and condition (g) that 
the transmission coefficient, r, is given by 

r=2(l+p), (13) 

which gives a simple relation between the reflection 



and transmission coefficient. Using the value of d 
found in (5), we have from (12) 

00 CO J 

S d m S _JL S rm sin qwr=kE (1 — p— g— pg) . (14) 

TO =0 r = l 2 7rr 

This scalar condition (together with (5)) will be used 
to evaluate p after (11) has been solved for the d n 's. 

4. Perturbation Method 

In order to use the perturbation method in the 
solution of the equations of the previous section, we 
must establish their explicit dependence upon k. 
Having done that, we shall solve the electrostatic 
problem (k=0) explicitly and use it as the unper- 
turbed case. 



Equation (11) depends upon k only in its term, 



k. 



lr 


■yjr 2 —k 2 


li- k - 


fen 


■\'n 2 —k 2 q 2 


n 1 &Y 



~W0-i)K?+--) 



_r / /k 2 k 2 q 2 \ 
n \ \r 2 n 2 J 

_!!_lZ_ /^ 2 . k 2 q 2 \ . 
n 2n \r 2 n 2 ) 

^jLaJ^l (1 g 2y> \ i 

~ n r 2 \nr WJ + 
Equation (14) depends upon k only in its term, l r . 

l —^2 :== - ir +2r~+ 

Neglecting all terms involving powers of k higher 
than the second, eq (11) can be written in matrix 
notation as follows: let 

-> 

d-=(d ,d ly d 2 , . . .) 

r=the matrix of elements T nr =- S Tn , n^O, 



=0, n=0, 



S— matrix of elements S, 

(as previously defined), 
C7=matrix of elements 



, _2 n 

rn 

K Jo 



n 



cos ry cos - ydy 
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=0 n or r=0. 

-> 
Now we add (5) to (11) and use the notation that 1 
is the vector 



1={1,0,0,0,...}. 

We assume that all summations go from to + o° . 



d=-{(T+k*U)S}(d)+E {l+p)l. 
If I is the identity matrix, we have 



Let 



{I+TS+VUS}(d)=E (l+p)l. 



T* m = S rm sin grr r=0,l,2, 

qw 



(15) 



(16) 



TJ* = — — T* r=l 2 

Then (14) becomes 

S *» S (T%+kV* Trrl )=kE (l-p-q-pq) (17a) 

w=0 r=0 

-> -> 

For shorthand, let us define vectors t* and u* whose 
components are: 



^ra jL—J ± rr> 
r=0 



..Hc^y-V TV* 



Then (17) can be rewritten as 



(t*+k 2 u*)-d=kE (l-p-q- P q). 



(17b) 
(17c) 

(18) 



Equations (16) and (18) form the system which 
we shall solve. First we shall solve (16) using the 

perturbation method to invert the matrix. Then 

-» 
we shall use the solution vector, d, in (18) to compute 
the reflection coefficient, p. 

The method will proceed as follows. We shall show 
that the electrostatic system 

{l+TS}(d)=R 

-» 
can be solved for arbitrary right-hand sides, R, in 



such a way that t* - d=0. When this is known, the 
full system can be solved up to terms of order k 2 by 
setting 



d=e+k 2 f. 

Equation (16) becomes 

{I+TS+k*US}&+k*f)=E (l+p)l. 
Thus, neglecting powers of k higher than the second, 

{I+TS}&+k^I+TS}(h+knUS}&) = E (l+p)l 

-> 
We now can find a unique vector, e, such that 

{i+rajoo=i 

-> 

(which implies that E (l-\-p)e solves the electro- 
-> -> -> 

static system) and 2*-e=0. Now, having found e, 
-» 

we compute {US} (e) and solve 

{I+TS}(f)=-{US}Ce) 

in such a way that t* -/=0. 

Having found d(=e+k 2 f)j we substitute this solu- 
tion vector into (18). 



E {l+ P ) 



Butt*-e=t*-j=0. 
Therefore, 



k 2 u* ■ e 



_ kE (l-p-q-pq) 
E (l+p) 

l — (q+ku*- e) 
l + (%+ku*-~e) 



(19) 



Our main result may now be stated as 

Theorem 1. For sufficiently small values of k, the 

reflection coefficient p is given by formula (19), where 

-> 
the vector u* has been defined by equation (17) and 

-> 
where the components oj the vector e are defined by 

-> 
equations (38) and (39); the vector e itself describes 
the solution of the electrostatic problem (k=0). 
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5. Electrostatic Case 
We have the system of linear equations 

{I+TS}(d)=B, 

and we wish to show that it can be solved for arbi- 
trary right-hand sides in such a way that the solution 

-> 
vector, d, satisfies the scalar condition 

To this end we first invert the matrix {I-\-TS}. 
The method we shall use is to produce a matrix, D, 
and a triangular matrix, A, such that 

{I+TS}{D}=A. 

Since A is triangular, its inverse can easily be 
computed; thus, the full solution is 

{I+TS}~ 1 =DA- 1 . 

First we shall construct the triangular matrix, A. 
To this end, let us consider the field in region I. It 
has the form 



El(x,y)=l£c n e Q cos-y 

n=0 <l 



(20) 



in the electrostatic case. If we relax the condition 
at + oo and allow the field to become exponentially 
large, new modes may be introduced. We shall 
refer to these as "unrealistic" modes. Let the iVtli 
unrealistic mode be denoted by E T N (x, y). It will 
increase exponentially for x^>™ , and the order 
of magnitude of the function at x^>co is deter- 
mined by the first term (n=—N) in the expansion: 



EMx,y)= S cSe **cos%. (21) 

n = —N 2 

We wish to determine E* N in such a manner that it 
satisfies all conditions for E y except for (7). This 
implies that, in region II, we have 



E^(x,y)=^2d^ n (x } y), 

71=0 

where, in the electrostatic case, 

<t>o(x ) y) = q+^ l S mo e mx cos my 

m=\ 

<t>n(x J y) = ^2S mn e mx cos my n^O 

m=l 

tn(x,y)=-^=JlmS mn e mx cos my 

(JX 771 = 1 



(22) 



all n. (23) 



Equations (21) and (22) provide us with an in- 
finitude of new fields. Solution of these will yield 
the matrices, D and A, which we seek. 

Matching (21) and (22) in the aperture x=0, 

\y\<z*, 

CO /TO CO 

X) c£cos-2/=2<tf»(0,2/). 

n = —N (£ n=0 

This implies that 



oK—JN 



Kl , 



(24) 



since 



and 



cos 



n / n\ 



cj> n (o,y)=cos-y, \y\<q7r. 



Let us define 



Is. 2 c o 



«=i ci\ 



Then 



, 2 L o 



ri>0 

n>N 
0<n<N 

n=0. 



dn=ci+c~ f for all n. 
We apply the derivative condition 

<)E l N {x J y) 1 ^~ N j* s m 



(25) 



dx 

*E%(x,y) 
dx 



Q_m = -N 

■■^]d^ n (x,y) 



S — (cm-Cm) cos ~ y=^d^ n (o, y) 
m=i q_ q 71=0 



TL 

Multiply by cos - y and integrate 



|y|<«pr- 

(26) 



n 2 (* Q7r °° n 

- (Cn—ei)=— S & m (o, y) cos - ydy. 

<L VJo m=o q 

Multiply this by ^ and add (25) 



00 2 r Qir r) 

2«=<fir+S « — Mo, y) cos - n ydy. (27) 

m=0 It to J (£ 



However, it is easily seen that (27) is merely 
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where 
and 



{I+TS}(d„)=2cZ 
d N =(f%,d?,...) 

c n = (2 c o> C-i> • • •> £-iV> 0, () . . .). 



(28) 



c N , will comprise the triangular matrix, A. 
Thus, we must do the following for each N: 



If for each N we can compute c N and d N , we will 
have found the matrices D and A. The solution 

vectors, d N , will comprise D; and the right-hand sides, 

fill comprise the triangular 
ms, we must do the followir 

(i) Compute the vector c#. 

(ii) Find the solution vector d N . 

Problem (i) may be solved as follows: We map 
the waveguide structure onto the strip \v\<7r in 
the w=u J riv plane. The Schwartz-Christ offel 
formula gives us the mapping im])licitly. We may 
state the result as 

Lemma 1. The mapping of the interior of the wave- 
guide in the z=(x-\-iy) plane (as drawn in fig. 1) 
onto the strip \v\ ^w of the w=(u J r%v) plane is given 

by 

\2 



Z = l0£ 



(vV^Hhi+V^+i) 2 



dw 
dz 



I e w +l 
\ q 2 e w +l 



(29) 



where the square roots are determined by 
arg (q 2 e WJ rl)^ = for w=iw J r2 log <z _1 + a- 



and cr^O, 



arg (q 2 e w -\-l)$ 



arg (e w + 1)1 = 
arg («*+l)*=| 



for w=iir-\-2 log <i lJ r<r 



and 0->O, 



for w^iir + T, r^O, 
for w=iT-\-r, r>0. 



The following points correspond to each other: 

w=iir—<n and z=iir— °° 

w=iir and z=iir 

w=iir+2 log q~ l and z=iirq 
w=iir-\-&> and 2 = i7rg + oo . 

Let f=6 % o-=e _T V Substituting into (29), 



f=o-- 



(V g 2 + a+Vl + cT) 



(Vg 2 +H-2Vi+<r) 2 (i-<z 2 ) 



(30) 



Q7T 



Figure 1. 



ay 



<£n = o 



a<£ n +k 2 <£ n = o 



• A. N 

I 9N=cos-y 



ay 



FlGUBE 2. 



Equation (30) shows that a can be expanded in a 
power series in f. The held E N can be written as 
the real part of the derivative of the potential, 



Let <K 



d 



e 



so that E N =He <$ N . We assert 



that $ v can be expanded in a scries of the type 

Nz 






e Q +c[%e 



(l-AO . 
(AT) Q 



+ ... +c' u ' 



(N) 






This form follows from (30) and from 

dz \ <j 2J r<r 

The residue theorem asserts that 






!) =ih$ tNf-Ut i^l^N. 



This reduces to the expression (31) for C { -} as an 
integral over a small circle around c=0, taken in the 
positive sense. 



niN) — 



N 



<z(i-<r) "<r=o 



(Vg 2 +cr+Vl + cr) 2 g 

(V?+^+<zVi+*) 2Z 



(31) 
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666696—63- 



This verifies that c { ^=0, />iV. 

Problem (ii) may be solved by using the fact that 

in the aperture each component of the vector d N is 
a coefficient of a Fourier series in the aperture. 



But 



m qrj. 



E N cos — ydy. 



(32) 



E N =TLe <t> w =Re -. Ne" 



, dw 
Tz 



Writing the full complex integral (and droppins,' the 
"Re" for now), 

1 f«*-'' dw f ~z --A 



. r w+loR G) eNW ( 



?« +e 



■ F(u>) 



} 



dw, 



where F(w) is the Schwartz-Chris toffel mapping 
given in (29). Using the substitution \=e w 



fi-Re 



N 



2tt qi 






J' 



X * 



(Vg 2 x+i+gVx+i) 
(V<z 2 x+i+Vx+T) 



d\+(l-q 2 ) 



(H)t 



J' 



A-i-f (V g 2 x+i+Vx+i) g 

(Vg 2 X+l+gV'x+l) 2ff ' 



d\ 



(33) 



The path of integration is around a circle of radius 

~2 in the positive sense, as illustrated in figure 3. 

It remains now to show that solutions of the 
system 



satisfy 



[I+TS}(d)=R 



X-plane 




Figure 3. 



The reason for this is that the components d m of d 

-> 
as well as the components of all the vectors d N are 
the Fourier coefficients of the derivative of a potential 

function Re ®—jr- (or, in general, of Re $ iV =^r- [sin 

Nv exp Nu\). In fact, the definition of the /* by 
(17b) shows that 



-r: 



tm(x,y)dy, 



where the \l/ m (x,y) are defined in the electrostatic 
case by (9a), (9b) with the additional condition 
k = 0. Now the d m are derived by the fact that 

2d m il' m (x'y) 

is the expansion for Re $ (or, more generally, Re $ N ) 
in region II. Therefore we have 



t*-d- 



~'2q 



-J. 



T d 2 v(o,y) d bv 



(0, ?tt) 



(34) 



dydx * dx (0t _ ffjr) 

But v is constant on both boundaries, so that there 
dv/dx = Q. Therefore, t*.d=0. 

6. Convergence 

It has been tacitly assumed in the previous 
computations that the convergence of the various 
infinite series is good enough to justify the operations 
performed. For example, the final formula for the 
reflection coefficient (19) involves an infinite series 

represented by the dot product, u*-e. It seems to 
be a rather difficult task to prove that our perturba- 
tion method leads to a convergent procedure for 
the computation of p, at least for sufficiently small k. 
All we shall do here is this: We shall prove that the 
approximation formula of Theorem 1 is meaningful. 
For this purpose, we must prove: 

Theorem 2: The infinite sum represented by the 
-> -> 
product u*-e converges absolutely. 

From equation (17) it can be seen that 



™_n 1 A (ITT I 



— Zqirr 



2 C q 

7T Jo 



in 7 

cos ry cos — ydy 



=— 2 S e, 



f ~ sin qirr C qir m i 1 

1 1 z] — ^— J o cos T v cos — ydy\ 

(35) 



We shall show first that the factor in brackets is 
less than a constant times 1/m for m large. This 
may be seen as follows: The factor in question is the 
mth Fourier coefficient of the function 



m^±{ s ^ r f +y) + sinr{i ;- y) y 
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The formula 

S 5 ^=-log(2cos(2//2)) 

r=l I 

shows that/(?/) is diflFerentiable excepl for isolated 
points in any finite interval. According to Whit- 
taker-Watson [1958], sec, 9.3, the mth Fourier 
coefficient of such a function has the ore lei- of mag- 
nitude of 1/m. 

In order to complete the proof of convergence of 
-» -> 

u* -e in (35), we must demonstrate the absolute 
convergence of 






(36) 



Each component of the vector e is a coefficient of 
the Fourier expansion of the field in the aperture. 



1 f«* 7- ™ 7 

Hy cos — my. 

9? J-q* 2 



(37) 



Inasmuch as E y vanishes for 2 7r< ^l?/l < C 7r > (37) can be 
written as 



cos — mlii 
<1 ' ' 



1 P F 

-=- — : -T- ( 6* +e « )dz. 

2irqiJ- 7ri dz \ J 



(38) 



We can convert this to a line integral in the ?/?-phmc 
using the Schwartz-Christoffel mapping given in (29). 






(V<z 2 ^+i+V^+i) 



2w 



+ 



1 



(!_gi)(l-f) 7 






— <iw y • 

(Vg^+i+gV^+i) 2 ™ f 

J 

(39) 



The path of integration, C, is illustrated in the 
following sketch (fig. 4): 

The first integral will converge if the path is 
altered in such a way that Ke (w)<(0, so we shall 
change the path into the one indicated by figure 5. 
Along the bottom path we may put 



w= 



-iri—s. 



0^s<< 



By rewriting the first integral in terms of s, and by 
using the inequalities 

(l-(/ 2 )^ (\-q 2 e- s )^+q(l-e- s ^ 

^(l-q 2 e-y + (l-e- s y*, 



we see easily that the resulting integral can be 
majorized by l/2irm. The same analysis goes 
through on the upper path of integration, and this 
si iows the absolute convergence of the first half 

of guv). 

For the second integral in (39) we can break up 
the path as indicated in figure 6. 

We see readily that the integrals along C 2 and C 3 
cancel each other out exactly. We are left, then, 
with integrals along C x and (7 4 . By putting 

W= iri-i-s, S S ^ log q' 2 

on c7 4 , and using the substitutions 

cos a=(l+e-*)i(l— q 2 )-*, 

cos/3=g(6'— 1)*(1— g 2 )~* (40) 



t=-\^ie-^)-Kl-e-T\ 



f=~2 «'(«'- 1)- } (1-<ZV)-*, 



(41) 
0gj8gjr/2 (42) 



»-r*s— J-' «^lCr") (43) 

we find that the second integral along (\ has the 
value 



( iy«+i riog a 2 

2ttk/ J 



(44) 



Applying a similar method to transform the second 
integral along C u we find, after combining the re- 
sults, the single integral 



-i) w+1 r h 

Kg Jo 






sin (2 my) ds, 



(45) 



w-plane 



Figure 4. 



log • 



^j. _ 



v-p!ane 



Figure 5. 
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w- plane 



Tri + log-^ 



-TTi+log-^ 



Figure 6. 

which is now our expression for the second integral 
in (39). It can be shown that (45) is majorized by 
a constant multiple of m~ 2/3 for large values of ra. 
The proof for this fact is somewhat tedious, and we 
shall give a sketch of it only. 

Because of (43), the integral in (45) is singular at 
s=log qr 2 , j8=7r/2, where ds/dy becomes infinite. 
To show that (45) exists and even becomes small for 
large m, we introduce the variable 



t=s— log q 2 

and show that t depends analytically on y in 
that (with a constant c) 



t=c 2 y 2/3J r higher terms in y 



1/3 



such 



(46) 



From (46), we can deduce that (45) behaves for 
large m like 

"* 5 sin 2 my 



f 

Jo 



T 



dy, 



(47) 



where 8 is a fixed, sufficiently small real number. 
In turn, (47) can be shown to be of the order of 
m~ 2/3 because of 

Combining all these results shows that u*-e con- 
verges absolutely, for we have shown that this sum 
behaves at worst like 



JZ 



1 



which is absolutely convergent. 



Figure 7. 

7. Concluding Remarks 

We have restricted our investigations to a first- 
order approximation of the value of the reflection 
coefficient. Higher-order approximations are avail- 
able, since we have the inverse of the matrix de- 
termining the electrostatic case. 

The method of solution used in this paper may be 
tedious, but it is sufficiently general to be applicable 
to a wider range of waveguide problems. For ex- 
ample, it can be used in finding the field in more com- 
plex waveguide structures such as those illustrated 
in the following sketches. 

In each of these problems there are three distinct 
regions, and this will give rise to two sets of reflection 
and transmission coefficients. (At first it might be 
thought that regions II and III in problem B could 
be combined into a single region, but elementary 
analysis shows that a field of this form is over- 
determined and can yield only the trivial solution.) 



The authors thank the referee for a large number 
of helpful suggestions. 
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